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ABSTRACT 

jn In a recently completed survey of the stellar populations properties of LINERS and 

LINER/HII Transition Objects (TOs), we have identified a numerous class of galactic 
nuclei which stand out because of their conspicuous 10 8-9 yr populations, traced by 
high order Balmer absorption lines and other stellar indices. These objects were called 
"Young-TOs" , since they all have TO-like emission line ratios. In this paper we extend 
, this previous work, which concentrated on the nuclear properties, by investigating 

the radial variations of spectral properties in Low Luminosity Active Galactic Nuclei 
(LL AGN) . Our analysis is based on high signal to noise long-slit spectra in the 3500- 
5500 A interval for a sample of 47 galaxies. The data probe distances of typically 
up to 850 pc from the nucleus with a resolution of ~ 100 pc (~ 1") and S/N ~ 
30. Stellar population gradients are mapped by the radial profiles of absorption line 
equivalent widths and continuum colours along the slit. These variations are further 
analyzed by means of a decomposition of each spectrum in terms of template galaxies 
representative of very young (< 10 7 yr), intermediate age (10 8-9 yr) and old (10 10 yr) 
stellar populations. 

rS I This study reveals that Young-TOs also differ from Old-TOs and Old-LINERs in 

terms of the spatial distributions of their stellar populations and dust. Specifically 
our main findings are: (1) Significant stellar population gradients are found almost ex- 
clusively in Young-TOs. (2) The intermediate age population of Young-TOs, although 
heavily concentrated in the nucleus, reaches distances of up to a few hundred pc from 
the nucleus. Nevertheless, the Half Width at Half Maximum of its brightness profile is 
more typically 100 pc or less. (3) Objects with predominantly old stellar populations 
present spatially homogeneous spectra, be they LINERs or TOs. (4) Young-TOs have 
much more dust in their central regions than other LLAGN. (5) The B-band luminosi- 
ties of the central < 1 Gyr population in Young-TOs are within an order of magnitude 

of Mb = —15, implying masses of order ~ 10 7 -10 8 M©. This population was 10-100 
times more luminous in its formation epoch, at which time young massive stars would 
have completely outshone any active nucleus, unless the AGN too was brighter in the 
past. 

Key words: galaxies: active - galaxies: Seyfert - galaxies: stellar content - galaxies: 
nuclei - galaxies: statistics 
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Low luminosity active galactic nuclei (LLAGN) are the most 
common form of activity in the nearby universe. Their prox- 
imity allows us to sample their properties on linear scales 
which are not accessible for farther AGN populations like 
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Seyferts and quasars. At optical-UV wavelengths, however, 
this advantage is compensated by the difficulty in isolating 
the light from these intrinsically weak nuclei out of a domi- 
nant stellar "background". 

In a series of papers, we have been working our own 
contribution to this field. In Paper I (Cid Fernandes etal. 
2004a) we analyzed ground based nuclear optical spec- 
tra of a sample of 51 LINERS and TOs, while Paper II 
(Gonzalez Delgado et al. 2004) complements this data set 
with archive HST/STIS spectra of 28 nearby LLAGN. These 
samples cover nearly half of the LLAGN in the survey of 
Ho, Filippenko & Sargent 1997 (hereafter HFS97). Our fo- 
cus throughout this series is on the stellar populations of 
LLAGN, with the ultimate goal of establishing the role of 
stellar processes in the physics of these objects. 

The main result of Papers I and II is that we have un- 
covered a very strong relation between the nuclear stellar 
population and the gas excitation, as measured by [OI]/Hq, 
the most important diagnostic line ratio in LLAGN. The 
relation is in the sense that virtually all systems containing 
strong populations of ~ 1 Gyr or less have [OI]/Ha < 0.25, 
while nuclei dominated by older stars span the full range in 
[OI]/Ha (up to nearly 1). In other words, virtually all sys- 
tems with relatively young stellar populations have TO-like 
emission line spectra, whereas older nuclei can have either 
TO or LINER-like line ratios. This finding lead us to intro- 
duce a combined stellar population and emission line clas- 
sification into 4 types: Young-TOs, Old-TOs, Old-LINERs 
and Young-LINERs. This latter class is extremely rare. 

This relation between line-ratios and stellar population 
is analogous to the one found in Seyfert 2s, where nuclei 
with strong circumnuclear starbursts tend to have relatively 
small values of [OIII]/H/3 and HeII/H/3, while systems dom- 
inated by old stars can reach larger values of these line ratios 
(Cid Fernandes etal. 2001 and references therein). It is thus 
tempting to interpret Young-TOs as low-luminosity analogs 
of starburst + Seyfert 2 composites, where the relatively 
low excitation is explained by the starburst contribution 
to H/3, which dilutes [OIII]/H/3 and HeII/H/3. Intriguingly, 
however, Young-TOs are substantially older than the star- 
bursts around Seyfert 2 nuclei, many of which are just a few 
Myr old, as deduced by the detection of O and WR stars. 
While Papers I and II revealed a surprisingly large number of 
systems containing 10 8 -10 9 yr populations, massive young 
stars of the type often found in Seyfert 2s seem to be rare in 
LLAGN. The analogy between Young-TOs and starburst + 
Seyfert 2 composites thus rests upon the hypothesis of the 
existence of a population of massive stars which remains es- 
sentially undetected at optical wavelengths. Clearly, further 
work is necessary to clarify the precise nature of the connec- 
tion between stellar and gaseous properties in LLAGN. 

One type of study which has been carried out for 
Seyferts is the mapping of stellar populations based on spa- 
tially resolved spectroscopy. Variations of absorption line 
equivalent widths (W\) and colours (C\) as a function of 
distance from the nucleus were mapped by means of long 
slit spectroscopy by Cid Fernandes, Storchi-Bergmann & 
Schmitt (1998), Boisson etal. (2000), Gonzalez Delgado, 
Heckman & Leitherer (2001), Joguet etal. (2001). These 
variations can be transformed into stellar population pro- 
files, as in the study by Raimann etal. (2003), who found 
that star-formation in starburst + Seyfert 2s composites, al- 



though concentrated in the central regions, is not confined 
to the nucleus, but spread over the inner ~ 1 kpc. Spatial 
gradients in spectral indices are also useful to detect the 
presence of a central continuum source, which dilutes the 
nuclear W\'s with respect to off- nuclear positions. Both a 
compact nuclear starburst and an AGN featureless contin- 
uum can produce this effect, but in Seyfert 2s the papers 
above showed that significant dilution only occurs when a 
starburst is present in the innermost extraction. 

Spatially resolved spectroscopy of LLAGN has so far 
been limited to relatively few studies (eg, Cid Fernan- 
des et al. 1998) . While these previous works advanced our 
comprehension of individual sources, the small number of 
objects, differences in spectral coverage, data quality and 
method of analysis prevents us from drawing general con- 
clusions about the radial distribution of stellar populations 
in LINERs and TOs. In this third paper we take advantage 
of our recently completed spectroscopic survey to extend 
this type of study to a large sample of LLAGN. Variations 
of spectral properties with distance from the nucleus are 
mapped with the general goal of investigating the relation 
between spatial gradients, emission line and nuclear stellar 
population properties. In particular, we aim at evaluating 
the spatial distribution of intermediate age populations, a 
distinguishing feature of Young-TOs. 

In g2j we describe the data set and present examples 
of our spatially resolved spectra. In we investigate the 
spatial variations of a set of spectral properties and quan- 
tify these gradients by means of suitable empirical indices. 
These gradients are further analyzed in ^with the goal of 
producing estimates of the sizes, luminosities, masses and 
extinction of the intermediate stellar population in the cen- 
tral regions of Young-TOs. These estimates provide useful 
hints on the past and future history of these sources. Finally, 
JO summarizes our results. 



2 DATA 

The data employed in this paper have been described in 
Paper I. Briefly, we have collected long-slit spectra in the 
3500-5500 A range for 60 galaxies selected out of the HFS97 
survey. Observations were carried out at the 2.5 m Nordic 
Optical Telescope with a 1" slit-width and the Kitt Peak 
National Observatory 2.1 m telescope with a 2" slit. Our 
survey differs from that of our mother sample in two main 
aspects: wavelength coverage and spatial resolution. The in- 
formation encoded in the region bluewards of 4200 A, not 
covered by HFS97, has been explored in previous papers in 
this series. Here we concentrate on the analysis of the spatial 
information in this data set. 

2.1 Extractions 

In order to map spectral gradients, spectra were extracted 
in several positions along the slit. Extractions for the KPNO 
spectra were made at every 2.34" (3 pixels) out to at least 
9 = ±4.7", but the seeing was 2-3" (FWHM). For the NOT 
spectra, which constitute 83% of the data analysed here, we 
have used 1.13" (6 pixels) long extractions out to at least 
6 = 4.5" from the nucleus in both directions. These narrow 
extractions approximately match the angular resolution of 
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our typical NOT observations, which were made under sub- 
arcsecond seeing. Outside this central region wider extrac- 
tions were used if necessary to ensure enough signal. 

The signal-to-noise ratio in each extraction was esti- 
mated from the rms fluctuation in the 4789-4839 A in- 
terval. Galaxies with (S/JV)a4soo ^ 15 at angular distances 
< 4.7" from the nucleus were deemed to have insufficient 
useful spatial coverage and discarded from the analysis. 
Our cleaned sample contains 47 objects, including 4 nor- 
mal galaxies and 1 Starburst nucleus. In the nuclear extrac- 
tions (S/N)\4soo varies between 31 and 88 with a median of 
51. Outside the nucleus, the median (S/N)\4,goo decreases 
from 45 at 9 = ±2.3" to 31 at ±4.5". The S/N in the 
4010-4060 A interval is typically 0.5(S/N) M8 oo- All 521 ex- 
tractions were dereddened by Galactic extinction using the 
Cardelli, Clayton & Mathis (1989) law and the Ab values 
of Schlegel, Finkbeiner & Davis (1998). We note that the 
Kitt Peak the observations (7 galaxies) were taken under 
non-photometric conditions. This, however affects only the 
absolute flux scale, not the shape of the spectrum, as we veri- 
fied comparing spectra of objects taken both in photometric 
and non-photometric nights. The single result reported in 
this paper which is affected by this problem is the luminos- 
ity of the central young population in NGC 404 f i|4.2.3ll . 
which is likely underestimated. 

The distances to the LLAGN in this sample vary be- 
tween d = 2.4 and 70.6 Mpc, with a median of 24.1 Mpc. 
At these distances, 9 = 4.5" corresponds to projected radii 
r = 52-1540 pc, with a median of 526 pc, while our nuclear 
extractions correspond to 11-204 pc in radius (median = 85 
pc). The spatial regions sampled by these observations are 
therefore smaller than the ones in our studies of Seyfert 2s 
(eg, Raimann etal. 2003), which sampled the inner few kpc 
with a resolution of ~ 300 pc. 



2.2 Sample properties 

Table [T] lists our sample, along with the useful spatial cov- 
erage in both angular (0 O ut) and linear (r out ) units, nuclear 
and off-nuclear S/N , linear scale, position angle and a sum- 
mary of spectral properties. 

The emission line classification from HFS97 is listed in 
column 8 of Table Q As in Papers I and II, we prefer to 
classify LLAGN as either strong or weak-[OI] emitters (col- 
umn 11), with a dividing line at [OI]/Hq = 0.25. These two 
classes differ only slightly from the LINER and TO classes of 
HFS97, and better represent the combined distributions of 
emission line and stellar population properties of LLAGN. 
Throughout this paper LINERs and TOs are used as syn- 
onyms of strong and weak-[OI] sources respectively. 

Paper I introduced a stellar population characterization 
scheme defined in terms of four classes: r\ = Y, I, I/O and 
O (column 9). The Y class denotes objects with a domi- 
nant young starburst. The only object in our sample which 
fits this class is the WR-galaxy NGC 3367, which is not 
a LLAGN but is kept in the analysis for comparison pur- 
poses. Nuclei with strong intermediate age (10 s -10 9 yr) pop- 
ulations, easily identified by High Order Balmer absorption 
Lines (HOBLs; H8A3889 and higher) and diluted metal lines, 
are classed as r\ — I, while nuclei dominated by old stars are 
attributed a r\ = O class, and rj — I/O denotes intermedi- 



ate cases. Not surprisingly, it is sometimes hard to decide 
where to fit a galaxy in this classification scheme. The best 
example of this sort of problem is NGC 772, which contains 
both young, intermediate age and old components (Paper 

I) . Despite the weak HOBLs in its spectrum, we chose to 
tag it as r; = I. 

A simpler (but still useful) classification scheme is to 
group rj = Y and I objects as "Young" and rj = I/O and O 
objects as "Old". As an objective criterion for this classifi- 
cation we use the value of the equivalent width of the Call 
K line in the nucleus: W^ uc < 15 A for Young systems and 
larger for Old ones (column 10). The use of this equivalent 
width as an indicator of the evolutionary status of the stel- 
lar population is justified because the AGN contribution to 
these continuum is these sources negligible (Papers I and 

II) . These two classes are paired with the [OI] /Ha class to 
produce our combined stellar population and emission line 
classification into Young/Old-TO/LINER, listed in the last 
column of Table Q 

Of the 42 LLAGN in our sample, 13 fit our definition 
of strong- [OI] sources and 29 are weak-[OI] sources, while 
the stellar populations types are split into 16 Young and 28 
Old systems. The combined emission line and stellar pop- 
ulation statistics are: 14 Young-TOs, 2 Young-LINERs, 11 
Old-LINERs and 15 Old-TOs. Note that Young-LINER is a 
practically non-existent category, as the overwhelming ma- 
jority of Young systems are weak-[OI] emitters. 

It is worth pointing out that Young-TOs in this sample 
are on-average closer than other LLAGN. The distances to 
Young-TOs span the d = 2.4-35.6 Mpc range, with a median 
of 16.8 Mpc, while for other LLAGN 14.3 < d < 70.6 Mpc, 
with a median of 31.6 Mpc. This tendency is already present 
in Paper I and in the HFS97 survey, from which we culled 
our sample. In principle one expects that radial variations of 
spectral properties due to the presence of a compact central 
source will be harder to detect for more distant objects, due 
to the increasing contribution of bulge light to the nuclear 
extraction. This potential difficulty, coupled with the trend 
discussed above may lead to a bias in the sense that radial 
gradients would be easier to detect in Young-TOs because 
of their smaller distances. We do not believe this effect has 
a strong impact on the conclusions of this paper, given that 
there is still a substantial overlap in distances of Young-TOs 
and other LLAGN. This issue is further discussed in «)3.1l and 

am 

2.3 Spatially resolved spectra: Examples and first 
impressions 

Figures and [5] illustrate spatially resolved spectra for a 
representative subset of the galaxies in our sample. Spatial 
gradients in spectral properties will be analyzed in detail in 
the remainder of this paper, but some results are evident 
from a simple visual inspection of these figures. 

(i) First, in objects like the Old-LINER NGC 315 the off- 
nuclear spectra look virtually identical to the nuclear spec- 
trum, implying a high spatial uniformity of the stellar popu- 
lations. The only noticeable gradient is in the emission lines, 
which are concentrated in the nucleus. 

(ii) Second, the strongest gradients are found in systems 
with conspicuous HOBLs (eg, NGC 4150, NGC 4569). As 
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Table 1. Col. (1): Galaxy name; Cols. (2) and (3): Useful angular and linear coverage. Col. (4): Angular scale. Cols. (5) and (6): S/N 
at 4800 A at nucleus an outer extractions. Col. (7). Slit position angle. Col. (8): Spectral type according to HFS97. Col. (9) Stellar 
population category (Paper I). Col. (10): Equivalent width of the Call K band at the nucleus, in A. Col. (11): W = Weak-[OI] (ie, 
[OI]/Hct < 0.25), S = Strong-[OI] ([OI]/Ha > 0.25). Col. (12): Combined emission line and stellar population class. Objects marked with 
a * were observed at KPNO. 



noted above, these are nearly all weak-[OI] sources. This 
combination of youngish stellar population and [OI]/Hcv < 
0.25 fits our definition of Young-TOs. 

(iii) Third, although HOBLs, when present, are stronger 
in the central extraction, they are not confined to the nu- 
cleus. This is clearly seen in the cases of NGC 4150 and NGC 
4569, where HOBLs still show up in extractions more than 



3" away from the nucleus. Given that the seeing in these 
observations was typically better than 1", we conclude that 
the "HOBLs region" is spatially extended. 

(iv) As is typical of LLAGN, emission lines are generally 
weak. In fact, many objects show no sign of important diag- 
nostic lines like H/3 and [OIII]A5007 even in the nucleus. The 
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Figure 1. Examples of spatially resolved spectra of LLAGN. Spectra have been normalized and shifted for clarity. The nuclear spectrum 
is drawn with a thicker line. Labels in the right indicate the angular distance from the nucleus. 



measurement of emission lines requires careful subtraction of 
the starlight, which we postpone to a future communication. 



3 STELLAR POPULATION GRADIENTS 

A convenient way to map spatio-spectral variations is 
to compute profiles of absorption features and contin- 
uum colours along the slit (eg, Cid Fernandes et al. 1998; 



Raimann et al. 2003). From Papers I and II we know that an 
AGN continuum contributes very little (if anything) to our 
ground based optical spectra. Any significant variation de- 
tected in these properties can thus be confidently attributed 
to variations in the stellar populations. 

In Paper I we have measured an extensive set of stel- 
lar population indices in different systems. In this paper we 
will use the equivalent widths of the Call K line (Wk), the 
G band (Wg), Mgl (Wm 9 ) and Wc (a "pseudo equivalent 
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Figure 2. As Fig.Q] 



width" centered in the continuum just to the blue of H9) 
plus the Cseeo = 3660/4020 and C5313 = 5313/4020 con- 
tinuum colours, all measured in Bica's system. The 4000 A 
break index of Balogh etal. (1999), D n (4000), is also used, 
but only for illustrative purposes. The Wc index works as a 
direct tracer of HOBLs: Spectra with clearly visible HOBLs 
all have Wc < 3.5 A, while spectra dominated by old popu- 
lations (~ 10 10 yr) have larger Wc due to a blend of metal 
lines. As shown in Paper I, Wk , which is a much stronger 



and thus more robust feature, is also a good (albeit indi- 
rect) tracer of the intermediate age populations responsible 
for the HOBLs. 



All these indices are highly correlated (Paper I). Their 
radial behaviors, however, need not be the same. For in- 
stance, a compact blue source such as young or intermediate 
age starburst should produce a larger dilution at the nucleus 
of the bluer indices, like Wk, than of the redder ones, such 
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Figure 3. Spatial variations of stellar population indices for NGC 
315. Top: Radial profiles of Wk (black, thick solid line), Wq (ma- 
genta, thin line), Wg (green, dotted line) and Wmq (red, dashed 
line). Note that Wc has been multiplied by 1.5 for clarity. Middle: 
Radial profile of the 3660/4020 (blue, dotted line) and 5313/4020 
(black, solid line) colours, and £> n (4000) (green, dashed line). The 
3660/4020 colour is multiplied by 2 in the plot. Bottom: Surface 
brightness at 4200 A (in flux units) along the slit, normalized to 
S(r = 0) = 1. The FWHM of the slit profile is marked as a thick 
line segment, and listed at the top right in both arcsec and pc. 
The dotted line shows the instrumental profile, corresponding to 
a star observed in the same night. Dotted and dashed vertical 
lines indicate projected distances of ±100 and 500 pc from the 
nucleus respectively. 



as Wm 9 . The comparison of the Wk and Wm s profiles may 
thus allow inferences about the nuclear stellar population. 

We have measured these indices automatically for all 
521 extractions analyzed in this work following the recipes 
outlined in Paper I for Bica's indices and Balogh et al. (1999) 
for D n (4000). This served as a further test of the objective 
pseudo continuum definition proposed in Paper I. After vi- 
sual inspection of the results, we have judged that only in 
18 spectra (3% of the total) the pseudo continuum deserved 
corrections. Uncertainties in all spectral indices were esti- 
mated by means of Monte Carlo simulations. The typical 
uncertainties at 9 = ±4.5", the outermost extractions in 
many of our sources, are 0.5-1 A for all equivalent widths, 
and 0.04 for C366O, C5313 and -D„(4000). Indices for the nu- 
clear extractions are 2-3 times more accurate. 



3.1 Radial Profiles of Stellar Indices 

Figures l3T ll2l show the variations of our seven stellar popula- 
tion indices with angular distance from the nucleus for some 
illustrative cases. The top panels show Wk (black, solid 
line), Wc (magenta, thin line), Wo (green, dotted line) and 
Wmq (red, dashed line). The middle panels show the C366O 
(blue, dotted line) and C5313 (black, solid line) colours, plus 
the D n (4000) profile (green, dashed line). The slit bright- 
ness profile S(r) at A = 4200 A is plotted in the bottom 
panel to give an idea of the light concentration. The thick 
line segment marks the FWHM of S(r); its value is listed 
in the top right in both angular and linear units. A stellar 
profile is also plotted to illustrate the spatial resolution. Ver- 
tical dotted and dashed lines indicate projected distances of 
±100 and ±500 pc from the nucleus respectively. 

The examples in fi gures l3t < 1 2 1 were chosen to illustrate 
the variety of radial profiles found in the sample. In a first 
cut, the W\ profiles may be grouped in three categories: 

(i) Flat (eg., NGC 305 and NGC 410), 

(ii) centrally peaked (eg., NGC 7742), 

(iii) "diluted" profiles (eg., NGC 3627, NGC 4569). 

Most objects studied here have either flat or diluted 
W\ profiles. In NGC 6951 and NGC 7742, the peaked ap- 
pearance of W\{r) is due to circum- nuclear star-forming 
rings which appear in our outermost extractions (Perez et al. 
2000). Outside these rings, the absorption lines rise up again, 
like in NGC 1097 and other ringed galaxies studied by Cid 
Fernandes etal. (1998). 

The main focus of our analysis throughout the rest of 
this paper will be on nature and properties of the source 
of dilution in LLAGN with diluted profiles. These profiles 
cannot be explained in terms of metallicity gradients, as this 
should produce peaked profiles. The drop in W\ towards the 
nucleus in these galaxies is thus clearly the result of dilution 
of the metallic features by a centrally concentrated stellar 
population which is younger than that a few arcseconds away 
from the nucleus. The most dramatic example of this effect 
is seen in the starburst galaxy NGC 3367, where the young 
starburst appears only in the three central extractions (fig- 
ure We note in passing that, at d — 43.6 Mpc, this 
galaxy is one of the most distant in our sample, well above 
the median distance of 27.9 Mpc. Yet, its W\ gradients are 
clearly mapped with our data, which shows that the worries 
raised in H2.2l about possible distance related biases and not 
justified in practice. Similar comments apply to NGC 5678 
(figure[71 d = 35.6 Mpc) and NGC 772 (figure[TTl d = 32.6 
Mpc). 

In LLAGN with diluted profiles, the diluting agent 
could in principle also be a young starburst, but, as shown 
in Papers I and II, in only ~ 10% of LLAGN such a young 
component contributes with more than 10% of the flux at 
4020 A in our ground-based nuclear spectra. For most ob- 
jects, the radial dilution is caused mainly by an intermedi- 
ate age population, which appears far more frequently and 
in much larger strengths. These populations are easily rec- 
ognized by their weak metal lines and deep HOBLs, as seen, 
for instance, in NGC 3627 and NGC 4569 (figures □ |U El 
and [Till. 

Figures H3Tll5l show the Wk profiles for all 47 galaxies in 
our sample, sorted in an increasing sequence of nuclear Wk 
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values. This ordering bears an excellent correspondence with 
the profile shapes: Of the first 19 galaxies (from NGC 3367 
to NGC 6500), at least 16 have diluted W K profiles. The ex- 
ceptions are NGC 2681, NGC 841 and possibly NGC 6500. 
From NGC 3166 onwards, ie, for W K UC > 15-16 A, profiles 
are either centrally peaked, or, more commonly, approxi- 
mately flat. This obvious link is examined in quantitative 
terms in the next section. 



3.2 Gradients in Equivalent Widths 

In order to quantify the spatial gradients seen in figures I3T 115I 
we define a radial dilution index 



Wf 



(i) 



which compares nuclear and mean off-nuclear equivalent 
widths. Flat W\ profiles should yield 8\ ~ 0, while S\ < 
correspond to centrally peaked profiles and Sx > to di- 
luted profiles. Furthermore, if the nuclear spectrum differs 
from that in off-nuclear extractions only by an extra con- 
tinuum source (or, more precisely, a source with negligible 
Wx), then S\ measures the fractional contribution of this 
source to the continuum at A (Cid Fernandes etal. 1998). 

Wf is defined as the average of W\{9) for extractions 
centered at \6\ between 2.2 and 4.7" from the nucleus. Note 
that for the NOT observations this definition excludes ex- 
tractions adjacent to the nucleus, which in some cases are 
contaminated by nuclear light due to seeing. The averaging 
is carried out weighting by the error in W\. The uncertain- 



ties in the dilution index were evaluated from standard er- 
ror propagation. Typical one sigma uncertainties in 8\ are 
0.1 for Wc and 0.04 for Wk, Wg and Wm 9 - We have also 
explored an alternative definition of W^ s in terms of extrac- 
tions between r = 250 and 750 pc from the nucleus, but this 
turned out to yield similar results, which further demon- 
strates that our conclusions are not significantly affected by 
potential distance-related biases ( H2.2I . 

Table [5] lists the resulting values of S\ . Gradients are 
considered to be significant whenever |<5k| > 10%, which 
corresponds to a ~ 2.5 sigma detection limit. According to 
this criterion, significantly diluted profiles (8k > 10%) occur 
in 13 of the 42 LLAGN in our sample, while only 3 have sig- 
nificantly peaked profiles (8k < —10%). Spatially homoge- 
neous stellar populations therefore prevail among LLAGN, 
accounting for ~ 60% of our sample. 

3.2.1 Relations between W\- gradients, emission line and 
nuclear stellar population properties 

In figure HBl we investigate the relation between dilution and 
nuclear stellar population by plotting 8\ against W^" for 
Wc, Wk, Wg and Wug- The vertical dotted lines in this 
plot are the same ones used in Paper I to approximately 
distinguish objects with significant intermediate age popula- 
tions (those with Wc < 3.5, W K < 15, W G < 9 and W Mg < 9 
A, which are classed as rj = /) from those dominated by 
older populations (n = I/O and O). Figure [TBI shows that 
these dividing lines also segregate objects with significant 
dilution from those without. Focusing on the W™ = 15 A 
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Figure 13. Radial profiles of the equivalent width of the Ca II K line for all galaxies in the sample. Dotted vertical lines mark distances 
of ±100 pc from the nucleus. A horizontal dashed line is drawn at Wk = 15 A for reference. The thick line segment in the bottom of 
each panel indicates the seeing, measured from the FWHM of star observed in the same night. Objects are sorted by the value of Wk 
at the nucleus (W K UC ), indicated in the bottom right corner of each panel. Galaxies in this figure have VFjJ- uc between 2.5 (bottom left 
panel) and 14.8 A (top right). 



limit, which separates Young from Old sources in our simple 
classification scheme, we find that 12 out of the 13 objects 
with 8k > 10% fall into the Young category, the exception 
being NGC 3245, which, with W K UC = 15.2 ±0.3 A, sits right 
at the border line between Young and Old sources. In other 
words, galaxies with significant radial gradients in their stel- 
lar populations contain intermediate age populations in their 
nuclei. The converse is also true, as at least 12 out of 16 



Young-LLAGN have diluted profiles. This is the same result 
found in figures I13H15I where we see that virtually every 
galaxy with W K UC < 15 A has a diluted W K profile. 

Since in Papers I and II we have shown that nearly 
all nuclei with weak metal absorption lines are weak- [01] 
emitters, we expect that the strong relation between S\ and 
W£ uc seen in figure Hoi translates to an equally strong rela- 
tion between 5\ and [OI]/Ha. This is confirmed in figure H7I 
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Figure 14. As Figure Il3l but for galaxies with Vl^ 110 between 14.9 and 17.8 A. 



which shows that all but one object with 8k > 10% have 
[01] /Hq < 0.25. Two other weak- [OI] nuclei, NGC 404 and 
NGC 4150, should probably be included in the list of sources 
with diluted profiles. NGC 404 is so close by (2.4 Mpc) that 
our outer useful extractions do not reach a probable rise in 
Wx's for larger radii, if this indeed happens in this dwarf 
galaxy. Dust effects may also be present, as indicated by the 
peak in the C5313 colour in the nucleus of NGC 404 (figure 
E2J). In NGC 4150 the strongest dilution is seen at 6 = +1.1" 
from the nucleus, and the rise in W\ seen in our last extrac- 
tions has a small weight in our definition of W£ s , resulting 
in a small 5\. This asymmetry is associated with the pro- 



nounced nuclear dust lane in this galaxy (Paper II), which 
is responsible for its asymmetric C5313 profile f figure 15^. 

The only strong- [OI] source with significant radial dilu- 
tion is NGC 5005 (6 K = 11 ± 3%, [OI]/Ha = 0.65). Given 
that this nucleus is classified as a LI. 9 by HFS97, it is con- 
ceivable that the dilution is caused by a nuclear featureless 
continuum, as found in spatially resolved spectroscopy of 
type 1 Seyferts (Cid Fernandes et al. 1998). However, none of 
the other 7 type 1 LLAGN in our sample exhibits significant 
dilution. Furthermore, HOBLs are clearly present in the nu- 
clear spectrum of NGC 5005, so we favor the interpretation 
that, as in other objects, dilution is caused mainly by a cen- 
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Figure 15. As Figure FL3l but for galaxies with W£ uc between 18 and 20 A. 



trally concentrated intermediate age population. As noted 
in Paper II, and confirmed by our radial dilution analysis, 
the contribution of a non-stellar continuum to our ground 
based spectra is negligible. Clear signatures of a featureless 
continuum in LLAGN are only found under the much higher 
spatial resolution of HST, and even then they are rare. 

We thus conclude that virtually all sources with radially 
diluted metal lines are weak-[OI] emitters. Note, however, 
that the converse is not true, as there are several weak- [01] 
objects with either flat or, more rarely, peaked W\ profiles. 
These non-diluted weak- [01] nuclei are dominated by old 



stellar populations, as deduced from their strong metal lines 
(figure [TBI Paper I). 

To summarize, combining the relations between dilu- 
tion, stellar population and emission line properties we find 
that significant stellar populations gradients are found al- 
most exclusively in Young-TOs, ie, objects with weak [01] 
and a conspicuous intermediate age nuclear stellar popula- 
tion. Old-TOs and Old-LINERs, on the other hand, tend to 
have spatially uniform stellar populations. These strong re- 
lations can be visualized comparing the location of different 
symbols in figures [TBI and [T71 
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Table 2. Cols. 2—5: Dilution of the equivalent widths Wc, Wk, Wg and Wug- Col. 6: Gradient in the C5313 colour (see equation |2*1. 
Cols. 7 and 8: HWHM of the flux profile along the slit, in angular and linear units. 



3.2.2 W\- gradients and the colour of the nuclear source 

Another result of the W\(r) analysis is that the spatial di- 
lution, when significant, tends to be larger for shorter wave- 
lengths, which implies that the diluting agent is bluer than 
the off-nuclear stellar population. This is illustrated in fig- 
ure 1181 where we plot the dilution in the K line (central 
A = 3930 A) against the dilution in W c (A = 3816 A), 
W G (A = 4301 A) and W Mg (A = 5176 A). For LLAGN 
with 8k ^ 10% the dilution follows a wavelength sequence: 
o~c > 8k > 8g > Sm 3 . (Deviations from this sequence are 
all within the uncertainties in S\). Some objects with clear 



gradients in K show little, if any, dilution in Mgl (eg, NGC 
3245 and NGC 6503). In NGC 772, NGC 4569 and other 
objects, the colour profiles confirm the existence of the blue 
nuclear component inferred from the behavior of Sx for dif- 
ferent lines. In others, however, C\(r) shows little variation 
(eg, NGC 3627) or even slightly redder colours in the nucleus 
(NGC 3245), contrary to the inference from the absorption 
line gradients. As discussed below, this apparent contradic- 
tion is due to dust in the central regions of these galaxies. 
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Figure 16. Radial gradients in four equivalent widths, measured from the comparison of nuclear and off-nuclear spectra. Different symbols 
correspond to Young-TOs (filled blue circles), Young-LINERs (open blue circles), Old-TOs (filled red triangles) and Old-LINERs (open 
red triangles). The star indicates the starburst galaxy NGC 3367. Crosses in the top right indicate mean error bars. Vertical dotted lines 
divide nuclei containing only old stars (large W\) from those with significant intermediate age populations (small W\). Note that NGC 
6951, which has very negative <5 A 's due to its star-forming ring, is outside all plot scales. 



3.3 Colour gradients and extinction 

Colour gradients carry information on the variations of stel- 
lar populations and extinction across a galaxy. Our C366O 
colour brackets the region containing the 4000 A break and 
Balmer jump, while C5313 is roughly equivalent to B — V. 
Because of the larger wavelength interval involved (5313 to 
4020 A) and the absence of spectral discontinuities in this 
range, C5313 is the more reddening sensitive of the two in- 



dices. One must nevertheless bear in mind that a Cs3i3(r) 
profile cannot be trivially transformed into an extinction 
profile without a simultaneous analysis of stellar population 
variations. 

Figures I19H21I show the C366O (dotted, blue line) and 
C5313 (solid, black line) colour profiles, also ordered accord- 
ing to W^. Centrally peaked Cs3i3(r) profiles are appar- 
ently rare among galaxies with W™ c < 15 A (figureEH , with 
exceptions (eg., NGC 404, NGC 718~and NGC 3245). This 
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Figure 19. Radial profiles of two continuum colours: C366O = 3660/4020 (dotted, blue line) and 5313/4020 (solid, black line). Note 
that the values of C5313 have been divided by 2 for plotting purposes. Dotted vertical lines mark distances of ±100 pc from the nucleus. 
Objects are sorted by the value of Wk at the nucleus, indicated in the bottom right corner of each panel. Galaxies in this figure have 
W™ c between 2.5 (bottom left panel) and 14.8 A (top right). 



type of profile appears more often in figures l^Ol and ETl which 
contain galaxies with W^- uc > 15 A. 

In order to examine colour gradients in more quan- 
titative terms we compute the ratio C^i 3 /C^ 13 between 
the values of C5313 in the nucleus and a mean off-nuclear 
colour, defined as the weighted average of extractions be- 
tween \0\ = 2.2 and 4.7" (as done for Wf in jOOJ. This 
ratio can be transformed into the index 



/ /-inuc \ 

S V = 5.98 log ^ (2 ) 

\°5313/ 

which measures by how many V-band magnitudes one has 
to deredden the nuclear spectrum to make it match the off- 
nuclear C5313 colour. The coefficient in this equation comes 
from assuming the Cardelli etal. (1989) extinction curve 
with Rv =3.1, which we do throughout this paper. 8v < 0, 
which indicates a bluening towards the nucleus, is hence- 
forth referred to as a "blue gradient" , while 5\r > is called 
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Figure 20. As Figure[HIl but for galaxies with W K UC between 14.9 and 17.8 A. 



a "red gradient" . Colour gradients were also examined by fit- 
ting the nuclear spectrum with a combination of off-nuclear 
spectra plus reddening, which yields the nuclear extinction 
relative to that of the off-nuclear extractions. This method 
gives essentially identical results to those based solely on the 
C5313 colour, with a mean off-set of just 0.04 mag and rms 
difference of 0.16 mag between the two Sv estimates. 

Figure HI compares S v with S K , W K UC and [OI]/Hq. 
The figure confirms that red gradients are more common for 
objects with W K UC > 15 A, ie, among Old-LLAGN, which 
also tend to have fiat W\ profiles. In fact, the plot shows 
that all objects with Wj^ 10 > 15 A, all strong- [OI] sources 



and all but one of the 8k < 10% objects have red gradients, 
the exception being NGC 4150, a Young-TO for which, as 
discussed above, 8k is underestimated. On the other hand, 
galaxies with significantly diluted W\ profiles, ~ 90% of 
which are Young-TOs, have both blue and red gradients. Of 
the 13 LLAGN with 8k > 10%, 8 have blue gradients and 5 
have red gradients, but note that in several of these objects 
the colour gradient is negligible, with \8v\ < 0.1 mag. 

Because colours per se do not disentangle intrinsic stel- 
lar population properties from extinction, these estimates 
of Sv can only be interpreted as actual spatial variations 
in dust content in the absence of stellar population varia- 
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Figure 21. As Figure Il9l but for galaxies with W K UC between 18 and 20 A. 



tions. This is a reasonable assumption for objects with rel- 
atively flat W\ profiles (and thus spatially uniform stellar 
populations) , which, as shown above, are essentially all Old- 
LINERs and Old-TOs. The red-gradients observed in these 
objects can thus be safely attributed to extinction gradi- 
ents. Note, however, that most objects with 5k ~ cluster 
around values of Sv of 0.1-0.4 mag, indicating that extinc- 
tion gradients tend to be small. 

The assumption of spatially uniform stellar populations 
breaks down for galaxies with diluted or peaked W\ profiles. 
In the latter case, one expects Sv to overestimate the extinc- 
tion gradient, as the reference off- nuclear extractions sam- 



ple a younger population than that present in the nucleus. 
This is clearly the case of NGC 6951, for which we obtain 
Sv = 1.2 mag, the largest value in the whole sample. This ef- 
fect is responsible for at least part of the trend of increasing 
of Sv as 5k becomes more negative (figure Conversely, 
when the nucleus contains a younger (and thus intrinsically 
bluer) population than off-nuclear positions, the resulting Sv 
should be regarded as a lower limit to the actual variation in 
Ay. From tl3.2l and figure IT51 we know that in galaxies with 
diluted W\ profiles the diluting source is intrinsically bluer 
than off-nuclear spectra, which should lead to blue gradi- 
ents. While some of these galaxies indeed have blue gradi- 
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ents, most (9/13) have negligible or slightly red gradients, 
which can only be understood in terms of a higher dust con- 
tent in the nucleus. Hence, contrary to the first impression 
derived from the relative rarity of centrally peaked C5313 
profiles among these sources, extinction gradients seem to 
be a common feature of Young-TOs. 

In summary, this empirical analysis shows that extinc- 
tion gradients are present in LLAGN of all kinds. In Old- 
LINERs and Old-TOs, which have spatially uniform stel- 
lar populations, these gradients are not huge, with 8v typ- 
ically smaller than 0.5 mag. Young-TOs, with their diluted 
W\ profiles, also have extinction gradients, but a quantita- 
tive assessment of their magnitude requires a more elaborate 
analysis, which we present in i|4.2.2l 



3.4 Slit Profiles 

The central intermediate age population which dilutes the 
equivalent widths of metal lines must cause an excess of flux 
with respect to the smooth surface brightness profile from 
the bulge of the host galaxy. Galaxies containing this extra 
central source should thus have sharper brightness profiles 
than those with more uniform stellar populations. 

In order to verify this prediction we have measured the 
Half Width at Half Maximum (HWHM) of the slit profiles, 
denoted by Rs ■ The results are listed in the last two columns 
of Table [5] in angular and linear scales, and graphically il- 
lustrated in figure 1251 The plot confirms that the most com- 
pact slit profiles occur among sources with diluted Wk pro- 
files. By extension of the relations between 5k, W k uc and 



[OI] /Ha, one expects these compact sources to be mostly 
Young-TOs, as confirmed in figures I^oT j and c. 

The slit profiles of Young-TOs suggests characteristic 
sizes of 50-100 pc for their central intermediate age popula- 
tion. This rough estimate suffers from two caveats. First, it 
is based on the total flux profile, which includes the bulge 
light. This issue is addressed in fQlbelow. Second, in several 
cases Rs corresponds to angular sizes of 1" or less (Table 
in which case seeing starts to dominate size estimates. 
In fact, the comparison of galaxy and stellar profiles in the 
bottom panels of figures l3T ll2l shows that while Old LLAGN 
have spatially resolved profiles, in Young TOs the inner S(r) 
profile is only marginally broader than the seeing disk, so 
Rs should be regarded as an upper limit for these objects. A 
more refined study of the inner morphology of LLAGN based 
on high resolution imaging is underway (Gonzalez Delgado 
etal., in preparation). 



4 ANALYSIS AND DISCUSSION 

Our spatially resolved spectra of LLAGN show that signif- 
icant stellar population gradients occur almost exclusively 
in Young-TOs. These gradients are caused mostly by an in- 
termediate age population (0.1-1 Gyr), although in a few 
cases a < 10 Myr nuclear starburst is also present (Papers 
I and II). The contribution of these stars to the total spec- 
trum increases towards the nucleus, causing the radial di- 
lution of metallic features. For consistence of notation, we 
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Figure 23. Half Width at Half Maximum of the flux distribution along the slit plotted against (a) the dilution in the K line, (b) the 
nuclear equivalent width of the K line, and (c) [OI]/Ha. Symbols as in figure Hoi 



NGC 4569 - Young-TO 



NGC 772 - Young-TO 



20 
15 

s 

js 10 



o 
o 
o 

<* 2 









: ; ; 








■ r t * 








1 1 1 


— i — i — i — 












1 — i — i — i — 




1 1 1 


I I ^-rrrT^T--' 


H h 

1 


1.3 " " 
107 pc ; 




-5 

distance from nucleus [arcsec] 



-5 5 

distance from nucleus [arcsec] 



Figure 10. As Figure^] but for NGC 4569. 



Figure 11. As Figure^] but for NGC 772. 
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Figure 12. As Figure |3 but for NGC 3367, a starburst galaxy 
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Figure 17. Radial dilution in the K line against the [OI]/Ha 
emission line ratio (extracted from HFS97). Symbols as in figure 
1161 The horizontal line shows the [OfJ/Ha = 0.25 taxonomical 
frontier which separates strong from weak-[OI] nuclei. Objects to 
the right of the vertical line at 8 k = 10% are those with significant 
dilution. 



hereafter denote this population the "Central Young Popu- 
lation" (CYP), where "young" means <1 Gyr-old. 

In this section we present estimates of the physical size, 
luminosity and extinction of the CYPs in Young-TOs. These 
estimates require separating the light from the CYP from 
that of older stars from the host's bulge, which in turn re- 
quires a more elaborated analysis than the eminently empiri- 
cal description of gradients presented in the previous section. 
Two methods were developed with this purpose. We close 
this section with a discussion on what these CYPs looked 
like in the past and what they might evolve to. 



4.1 Fits of the Equivalent Width profiles 

A rough estimate of the size of the region responsible for the 
radial dilution of metal lines may be obtained by evaluating 
at which distance from the nucleus WK{r) crosses the divid- 
ing line at Wk = 15 A, which characterizes the transition 
from "Young" to "Old" stellar populations in our simple 
classification scheme. This is not always possible, either be- 
cause Wk sometimes does not raise above this threshold in 
the whole region analyzed (eg., NGC 4569, figure HUl or be- 
cause of asymmetries or oscillations in the Wk profile (eg., 
NGC 5678, figure0. For the objects where this analysis was 
possible, we estimate radii between ~ 100 and 300 pc. 

A more formal estimate may be obtained fitting the W\ 
profiles. A two-components model was build for this pur- 
pose. We assume that W\ (r) results from the superposition 
of a "background" component with a W\(r) = Wa(oo) flat 



K 40 




Figure 18. Dilution in Wc (diamonds), Wq (crosses) and Wmq 
(filled circles) against the dilution in Wk- The diagonal line sim- 
ply marks y = x. All sources in this plot are LLAGN. Error bars 
have been omitted for clarity. 
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profile and a diluting component with negligible W\, and 
whose fractional contribution f(r) to the total continuum 
at wavelength A and position r follows a bell-shape radial 
distribution. The resulting model is expressed by 



Wx(r)=W x (oo) [l-/(r)] = 



oo) 



1 - 



1 + {r/aw) 2 



(3) 



where A A = [Wa(oo) - Wx(0)]/Wx(oo) and Wx(oo) are the 
analytical equivalents of Sx and respectively (see equa- 
tion 0. aw is the HWHM of the f(r) profile, a size scale 
which should not be confused with the HWHM of surface 
brightness profile associated with the central diluting com- 
ponent. The latter quantity, which we denote by Rw, must 
be evaluated from the product of S(r) and f(r). 

We have fitted this model to the Wk profiles of 15 
LLAGN: NGC 404, NGC 4150 plus the 13 LLAGN with 
5k > 10% (ie, those with significant dilution). The results 
are reported in Table |3] The fits are generally good, as illus- 
trated in figure |2"11 The dilution factors obtained from the 
fits are larger than the ones measured through equation Q, 
with A k ~ 1.35k typically. This happens because in most 
cases our operational definition of Wjc includes part of the 
rising portion of the Wk(t) curve, while equation © fits an 
asymptotic value. Interestingly, we find Ak = 23 and 29% 
for NGC 404 and NGC 4150 respectively, two Young-TOs for 
which 5k fails to detect significant dilution but, according 
to the qualitative considerations in H3.2I should be included 
in the list of sources with diluted Wk profiles. 

The values of aw range from ~ 30 to 400 pc, with a 
median of 172 pc, in agreement with the cruder estimates 
based on the size of the W K < 15 A region. These val- 
ues are larger than the HWHM of S(r), which spans the 
Rs — 20-170 pc range, with a median of 93 pc for this sub- 
set of galaxies (Table |5J. In other words, f(r) is broader 
than S(r). Therefore, in practice the HWHM of the light 
distribution associated with the CYP is dictated more by 
the slit profile than by the f(r) deduced from the WK{r) 
fits. The S x / profiles yield Rw = 17 to 126 pc (median 
= 67 pc), just slightly smaller than Rs- Hence, although it 
is clear that these CYPs often extend to more than 100 pc 
from the centre (as demonstrated by the detection of HOBLs 
well outside the nucleus; eg., figures Q and |2J , most of their 
light is concentrated within r < 100 pc. The relatively little 
light from the outer (r > Rw) parts of this distribution is 
enough to compete with the flux from the host's bulge, pro- 
ducing diluted Wx profiles on scales aw substantially larger 
than Rw- 

It is important to emphasize that in angular units, the 
median Rw corresponds to just 0.8". Hence, although we 
are able to resolve the wings of the light profile of CYPs, 
seeing prevents us from adequately sampling their core. Our 
estimates of Rw should thus be regarded as upper limits to 
the actual CYP radius. Indeed, high resolution images of 
a few Young-TOs reveal structures on scales smaller than 
the ones we are able to trace with our ~ 1" resolution. For 
example, NGC 4569 is known to have a very strong and com- 
pact nuclear source. Maoz etal. (1996) found, based on an 
HST/FOC image at 2300 A, that the emission of this galaxy 
is composed of a bright unresolved nuclear point source 
and some faint extended emission 0.65" south of the nu- 
cleus. Similar observations by Barth etal. (1998), done with 



HST/WFPC2 at 2200 A, find that the nucleus is slightly 
resolved along PA = 20°, with a dimension of 0.16" x 0.11". 



4.2 Template decomposition 

4.2.1 Method 

An alternative and more complete way to analyze gradi- 
ents in stellar populations is to model each extraction in 
terms of a superposition of spectra of well understood stellar 
populations. This can be achieved by means of the empir- 
ical starlight modeling scheme introduced in Paper I. The 
method consists of fitting a given spectrum with a combina- 
tion of five non-active galaxies from our comparison sample, 
whose spectra represent stellar population classes r\ = Y 
(NGC 3367), / (NGC 205), I/O and O (NGC 221, NGC 
1023 and NGC 2950). The code outputs the fractional con- 
tribution of these components to the flux at 4020 A, ex- 
pressed as a population vector x = (xy,xi,Xo), where the 
T) = I/O and O components are grouped in xo for concise- 
ness. The code also fits the extinction Av, modeled as due 
to an uniform dust screen with Av up to 4 mag. Regions 
around emission lines are masked out in the comparison of 
model and observed spectra. Paper I shows that this method 
provides excellent fits to the spectra. Unlike in Papers I and 
II, we have dereddened the template galaxies by their intrin- 
sic extinction derived by method described in Cid Fernandes 
etal. (2004b). Only NGC 3367 and NGC 205 are found to 
have significant extinction, both with Av ~ 0.9 mag. These 
corrections were applied because of our interest in estimat- 
ing the extinction and its radial variations in LLAGN. 

We have applied this method to all nuclear and off- 
nuclear spectra analyzed in this paper, thereby producing 
stellar population and extinction profiles. The spectral fits 
are of similar quality to those exemplified in Paper I. The 
median fractional difference between model and observed 
spectra for all extractions is 4.5%, which is acceptable con- 
sidering a median noise-to-signal ratio of 6% at 4000 A and 
3% at 4800 A. 

Examples of the resulting x(r) and Ay(r) are illustrated 
in figure 1251 The population vector in these plots is grouped 
into a predominantly old component, xo, (dotted red line) 
and a young + intermediate age component, xy+i = xy+xi 
(solid blue line), representing the combined strengths of the 
NGC 3367 and NGC 205-like components. This coarse 2- 
components description of the stellar population matches 
our Young/Old classification scheme. The xy+i{r) fraction, 
in particular, is used to map the CYP. We further plot xy (r) 
as a dashed line to illustrate that xy+i is actually dominated 
by the intermediate age population. As found in Papers I 
and II, young starbursts are generally weak or absent in 
LLAGN, although off- nuclear star-formation occurs in a few 
cases, as in NGC 6951. This Old-TO provides a good exam- 
ple of the power of the method (figure 1231 : L Its well known 
star-forming ring (Perez et al. 2000) , at r ~ 4" ~ 500 pc, is 
nicely mapped by the xy+i profile and its associated bright- 
ness distribution (bottom panel), obtained from the multi- 
plication of xy+i(r) by the slit profile S(r). Notice also the 
rise in extinction in the ring, the presence of an intermediate 
age component throughout the observed region, particularly 
in the ring, and the prevalence of an old, bulge-like popu- 
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Figure 24. Top: Examples of the Wk{t) fits. Bottom: Crosses show the normalized slit brightness profile. The solid blue line shows the 
brightness profile of the diluting component inferred from the Wk(t) fits. Labels indicate the HWHM of the total brightness profile (Rs) 
and the HWHM of the diluting source (Rw), also indicated by the thick horizontal line-segment. Vertical dotted lines mark projected 
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lation in the nucleus, which accounts for over 70% of the 
light. 

Figure l2"6l shows the x(r) profiles for all 42 LLAGN in 
our sample. As in previous plots, galaxies are ordered from 
bottom-right to top- left in an increasing sequence of W K UC . 
The plot confirms that the spectral gradients identified in 
S0are indeed associated with a centrally concentrated inter- 
mediate age stellar population, plus, in a few cases, a young 
starburst (eg., NGC NGC 772). This can be seen by the 
peaked x Y +i profiles from NGC 4569 to NGC 6500 in figure 
1261 Conversely, the old stellar component, mapped by xo(r) 
in these plots, bears a clear similarity in shape with the 
W\ profiles: Galaxies with diluted lines have diluted xo(r) 
profiles. Similarly, the spatial homogeneity of stellar popula- 
tions inferred from the flat W\(r) in galaxies like NGC 266 
and most others in the right half of figure 12^1 is confirmed by 
equally flat xo profiles, while peaked W\ profiles map onto 
peaked xo profiles (eg, NGC 1161). 

Hence, to first order, the x(r) profiles obtained from the 



template decomposition merely map the W\ variations onto 
the stellar population space spanned by our normal galaxy 
base. In fact, this relation is so strong that the equation 



xo = (0.068 ± 0.001) [A] - (0.35 ± 0.02) (4) 



transforms Wk into xo to within better than 0.1 rms for 
all 521 spectra. Plugging our W^ uc = 15 A dividing line in 
this equation we find that the transition from Young to Old 
stellar population occurs around xo ~ 2/3, or, equivalently, 
xy+i ~ 1/3. We thus conclude that CYPs which account for 
< 1/3 of the optical light would not be recognized as such 
in our data. Indeed, of the 15 LLAGN with CYPs detected 
through the radial dilution of Wk only 2 have xy+i < 1/3: 
NGC 4826 (W K UC = 14.4 A and x Y +i = 0.25) and NGC 
3245 (W K UC = 15.2 A and x Y +i = 0.23). 
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Results of the Wk (r) fits 



NGC W K {oo) [A] A K [%] a w ["] a w [pc] ij w ["] R w [pc] 
(1) (2) (3) (4) (5) (6) (7) 



0404 


12.2 


± 


0.4 


23 


± 


4 


3.7 


± 


0.5 


43 ± 6 


1. 


.7 


± 


0. 


1 


20 ± 1 


0718 


17.9 


± 


1. 1 


28 


± 


6 


3.6 


± 


0.6 


374 ± 65 


1. 


.2 


± 


0. 


1 


126 ± 5 


0772 


19.9 


± 


0.3 


42 


± 


3 


1.1 


± 


0.1 


172 ± 20 





.5 


± 


0. 


1 


86 ± 3 


3245 


18.2 


± 


0.3 


17 


± 


2 


1.0 


± 


0.3 


107 ± 29 





.6 


± 


0. 


1 


67 ± 8 


3627 


15.5 


± 


0.3 


25 


± 


3 


1.4 


± 


0.3 


44 ± 9 





.8 


± 


0. 


1 


25 ± 2 


3705 


20.4 


± 


0.7 


28 


± 


4 


1.7 


± 


0.5 


139 ± 39 





.6 


± 


0. 


1 


52 ± 3 


4150 


15.7 


± 


1.3 


29 


± 


7 


4.2 


± 


0.7 


199 ± 32 


1. 


.8 


± 


0. 


1 


82 ± 3 


4569 


22.6 


± 


2.4 


78 


± 


2 


4.5 


± 


0.5 


365 ± 37 





.6 


± 


0. 


1 


52 ± 3 


4736 


15.6 


± 


0.2 


17 


± 


2 


2.3 


± 


0.6 


47 ± 12 





.9 


± 


0. 


1 


19 ± 1 


4826 


17.0 


± 


0.6 


15 


± 


4 


2.0 


± 


0.7 


39 ± 14 





.8 


± 


0. 


1 


17 ± 2 


5005 


16.5 


± 


0.3 


15 


± 


3 


1.7 


± 


0.2 


172 ± 24 





.8 


± 


0. 


1 


85 ± 4 


5377 


18.6 


± 


1.0 


52 


± 


2 


2.6 


± 


0.4 


388 ± 59 





.7 


± 


0. 


1 


106 ± 2 


5678 


13.2 


± 


0.7 


33 


± 


4 


1.1 


± 


0.7 


190 ± 115 


0. 


.6 


± 


0. 


2 


110 ± 34 


5921 


18.2 


± 


1.7 


38 


± 


6 


2.0 


± 


0.5 


244 ± 63 





.8 


± 


0. 


1 


96 ± 7 


6503 


14.4 


± 


0.6 


34 


± 


4 


1.1 


± 


0.1 


31 ± 13 





.6 


± 


0. 


1 


17 ± 3 



Table 3. Results of the W K (r) fits for LLAGN with diluted W K profiles. 



4-2.2 Extinction profiles 

Our empirical analysis of colour and equivalent width gradi- 
ents in M. V indicates that extinction gradients are generally 
small in Old-LLAGN, while for Young systems we could only 
reach the qualitative conclusion that extinction variations 
must occur. A much more refined analysis is possible with 
the template decomposition method, which produces quan- 
titative estimates of both gradients and absolute values of 
the extinction. 

The Ay(r) profiles derived by this method are presented 
in figure 1771 for our 42 LLAGN. The first result which strikes 
the eye in this plot is the obvious asymmetry between galax- 
ies in the left and right halves of the figure, which, given 
the ordering according to WJ UC , essentially correspond to 
Young and Old systems respectively. The extinction pro- 
files of Young- LLAGN are substantially more complex than 
those of Old-LLAGN, which are often approximately flat. In 
both cases, extinction gradients, when present, are generally 
in the sense of producing centrally peaked Ay profiles, indi- 
cating a higher concentration of dust in the central regions. 
It is nevertheless clear that other types of dust distribution 
exist, as in NGC 4150 and NGC 4826, whose asymmetric 
Ay(r) curves indicate the presence of off- nuclear dust-lanes. 

A second and even more obvious result from figure 1771 
is that there is a clear offset in the absolute values of Ay 
between Young and Old systems. The statistics of Ay reflect 
this difference. Averaging Ay (r) over all extractions for each 
galaxy, we obtain a median spatially-averaged extinction of 
0.42 for our 16 Young-LLAGN, compared to 0.11 for the 
26 Old-LLAGN. A similar off-set is found considering only 
the nuclear extractions, which have median Ay(0) = 0.62 
and 0.21, respectively. Young-LLAGN, ~ 90% of which 
are Young-TOs, are therefore ~ 3 times dustier than Old- 
LLAGN. The clearest exception to this strong correlation is 
NGC 4438. The high concentration of dust inferred from Ay 
profile of this Old-LINER is associated with the pronounced 
nuclear dust lane seen in HST images (Kenney & Yale 2002) . 

The Balmer decrement measurements of HFS97 lend 
further support to interpretation that Young-TOs have a 
higher dust content than other LLAGN. Using their tab- 



ulated values for objects in our sample, we find a median 
Ha/H/3 of 4.6 for Young-TOs and 3.1 for other LLAGN. 
We can extend this analysis to the whole HFS97 sample 
using their measurements of the G-band equivalent width 
and classifying LLAGN into Young or Old adopting a W(G- 
band) = 4 A dividing line, which is roughly equivalent to our 
Young/Old division at W K = 15 A (Paper I). The 27 Young- 
TOs in this larger sample have a median Ha/H/3 = 4.5, while 
for the other 116 LLAGN this ratio is 3.2. 

We thus conclude that all evidence points towards a 
scenario where Young-TOs are the dustier members of the 
LLAGN family. 



4-2.3 Sizes and luminosities of the CYPs 

The population vector derived through the template decom- 
position analysis may be combined with the slit-profiles to 
produce one- dimensional surface-brightness profiles of the 
different stellar populations in our galaxies, as illustrated 
in the bottom panels of figure 1251 In what follows we use 
this method to estimate sizes and luminosities of the CYPs, 
represented by the Scyp(t) = S(r) x xy+i{r) profile. This 
method differs from the one in t|4. II in two aspects: (1) in- 
stead of assuming a functional form for the light fraction 
associated to the CYP we derive this fraction empirically 
from the template decomposition; and (2) all the spectrum 
is used, as opposed to a single equivalent width. 

Figure 1281 shows the total slit profile S(r) (thin black 
line) , and its decomposition into Young (thick blue line) and 
Old (dotted red) components for our 42 LLAGN. The plot 
shows that the young components dominate the light in the 
inner ~ 100 pc from NGC 4569 up to NGC 718 {W 1 ^ = 
13.1 A), becoming fainter than the inner old population as 
VK]J- UC increases, until it eventually "vanishes" from NGC 
7177 onwards (W^ uc > 16.6 A). Note that, unlike all other 
profiles in this paper, figure l2^1 uses a linear scale for r, which 
emphasizes the compactness of the CYPs in Young-TOs. 

We estimate the radius of the CYPs from the HWHM 
of the ScYp(r) profiles. Table 2] presents our results. As for 
the Wk(t) fits, we obtain xy+i profiles which are broader 
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(a) NGC 5678 (b) NGC 3627 (c) NGC 6951 
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Figure 25. Examples of the results of the spectral fit of the spatially resolved spectra of four LLAGN using a base of template galaxies. 
Top: Radial profile of the population vector. Old stellar populations (xj/q + xq) are represented by dotted (red) lines, young plus 
intermediate and age populations (xy+i = xy by a solid (blue) line and young starbursts (xy) by a dashed line. Middle: Radial 

profile of the extinction Ay(r). Bottom: Crosses show the observed surface brightness profile along the slit, normalized to its peak value. 
The solid (blue) and dotted (red) lines represent the profiles associated with the young + intermediate and old components respectively. 
Vertical dotted lines indicate projected distances of ±100 pc from the nucleus. 



than S(r), so Rcyp is close to Rs (Table |5J. The values of 
Rcyp are in good agreement with Rw (Table [3J, which is 
the equivalent CYP radius in the Wk(t) fits. Again, these 
estimates should be regarded as upper limits given that the 
angular sizes are limited by our spatial resolution. 

The luminosity associated with the CYPs was estimated 



integrating ScYp(r) within \r\ < 5Rs- The integration is 
performed in half-rings of area nrdr, ie., extrapolating our 
ID profiles to 2D. Table |U lists both the total and CYP 
luminosities. Numbers in between parentheses correspond 
to luminosities corrected by intrinsic extinction using the 
modeled Av profiles. The resulting dereddened CYP lumi- 
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Figure 26. Results of the template decomposition for all 42 LLAGN. Plots are ordered according to the value of W™ c , from small 
values in the bottom-left to large values in the top right. Dotted, red lines correspond to xq\ solid blue lines correspond to xy+i and 
dashed lines to xy 



nosities at 4020 A range from Lcyp ~ 10 3 - 3 to 10 5 5 L Q /A, 
with a median of 10 4 ' 3 Lq/A. Expressed in more conven- 
tional units, this roughly corresponds to a range in B-band 
absolute magnitudes 1 from ~ —12.2 to -17.7, with a median 
Mb = —14.7. Given the uncertainties in absolute flux cali- 



bration, extinction correction and extrapolation from ID to 
2D profiles, these values should be taken as order of magni- 
tude estimates. Yet, they are precise enough for the general 
considerations we present next. 



1 We use a Mg — 2.51ogL4Q20 — 3.96 conversion, for L4020 m 
units of Ij0 /A, derived from the Starburst99 models. 
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Figure 27. Extinction profiles obtained from the template decomposition analysis. The ordering of the galaxies is as in Fig. I2ril 



4.3 Discussion: The past and future of Young-TOs 

Naturally, the intermediate age stars which typify the CYPs 
of Young-TOs have been younger in the past and will get 
older in the future. Their current age and luminosity can 
be used, with the aid of evolutionary synthesis models, to 
predict what these objects looked like in their early days and 
what they will eventually become. 

For simplicity, lets assume that CYPs formed in instan- 
taneous bursts 10 8 -10 9 yr ago. From the Starburst99 models 
of Leitherer etal. (1999) one infers that these CYPs were ~ 



10 to 100 times more luminous in the optical in their first 
Myrs of life. Since the old stellar population has not changed 
substantially over this period, the CYPs would be much eas- 
ier to detect back then. The weakest CYPs recognized as 
such in our sample (ie, those with W™ c < 15 A) presently 
account for xy+i ~ 33% of the nuclear light at A4020. Scal- 
ing their present luminosity by factors of 10-100 would raise 
this fraction to 83-98%, which shows that they would com- 
pletely outshine the bulge light, and the optical continuum 
would be essentially identical to that of a starburst galaxy. 



Stellar Populations of Low Luminosity Active Galactic 27 




-0.5 0.5 -0.5 0.5 -0.5 0.5 -0.5 0.5 -0.5 0.5 -0.5 0.5 



distance from nucleus [kpc] 



Figure 28. Normalized slit brightness profiles (thin black line), decomposed into young (thick solid blue line) and old components 
(dotted red). The ordering of the galaxies is as in Fig. 1261 



Recall, however, that Young-TOs are dusty, so these lumi- 
nous infant CYPs could be substantially reddened and thus 
powerful far-IR sources, particularly if they had even more 
dust (and gas) in their early phases. 

The hot, massive stars in these early phase would have 
a large impact in the ionizing photon field. The Ha to A4020 
flux ratio for young starbursts is of order 1000 A (Leitherer 
etal. 1999). Currently, CYPs have L XA020 ~ 10 4 3 L A" 1 
(Table 2J, which scaled back to t = yields Ha luminosi- 



ties of order 10 42 ergs -1 , more than two orders of magni- 
tude larger than those currently observed in Young-TOs and 
LLAGN in general, which range from 10 38 to 10 40 ergs" 1 
(HFS97). In terms of Lhc, they would rank among pow- 
erful starburst nuclei and Seyferts. Clearly, these objects 
would definitely not be classified as "Low Luminosity" in 
their youth. It is not clear whether they would be classified 
as AGN either! Unless the AGN too was much brighter in 
the past, these objects would surely look like starbursts 
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CYP Sizes & Luminosities from Template Decomposition 



NGC 


Rcyp ["] 


Rcyp [pc] 


lo 


gL tot 


log L C YP 


(1) 


(2) 


(3) 




(4) 


(5) 


0404* 


2.5 


29 


3.77 


(4.08) 


3.54 (3.87) 


0718 


1.3 


133 


5.14 


(5.17) 


4.59 (4.64) 


0772 


0.5 


79 


5.06 


(5.21) 


4.10 (4.32) 


3245 


0.7 


74 


5.15 


(5.19) 


4.18 (4.22) 


3367 


0.4 


90 


5.00 


(5.35) 


4.91 (5.28) 


3627 


0.9 


30 


4.31 


(4.60) 


3.93 (4.26) 


3705 


0.6 


48 


4.24 


(4.44) 


3.48 (3.72) 


4150 


2.1 


98 


4.53 


(5.00) 


4.15 (4.70) 


4569 


0.6 


53 


5.20 


(5.53) 


5.19 (5.50) 


4736 


1.0 


21 


4.64 


(4.77) 


4.20 (4.34) 


4826 


0.9 


19 


3.89 


(4.05) 


3.11 (3.36) 


5005 


0.9 


98 


5.15 


(5.45) 


4.55 (4.93) 


5377 


0.7 


104 


5.11 


(5.36) 


4.85 (5.10) 


5678 


0.8 


135 


4.89 


(5.53) 


4.60 (5.27) 


5921 


0.8 


95 


5.03 


(5.21) 


4.63 (4.83) 


6503 


0.7 


19 


3.30 


(3.61) 


3.01 (3.33) 



Table 4. CYP size and luminosities estimates from the template decomposition analysis. Columns 4 and 5 give the total and CYP 
monochromatic luminosities at 4020 A integrated along the slit and extrapolated to 2D, in units of LqA -1 . Numbers in between 
parentheses are the dereddened luminosities. * = Observed under non-photometric conditions. 



Simple stellar populations of ages between t ~ 10 s and 
10 9 yr have mass-to-light ratios at A4020 of ~ 500 to 5000 
Mq Lq 1 A in the solar metallicity models of Bruzual & Char- 
lot (2003). For a median CYP luminosity of 10 4 3 L Q A" 1 
(Table El, this implies CYP masses M C yp ~ 10 7 -10 8 Mq. 
Star-formation has either ceased long ago or proceeds at 
a residual level in CYPs, otherwise they would look much 
younger. It is thus reasonable to suppose that these stars 
formed over a period of time whose length is a fraction of 
their current age. For star-formation time scales of 10 7 -10 8 
yr, the typical star- formation rate was of order 1 Mq yr" 1 . 
These are clearly very rough estimates, but they serve to set 
the scale of the CYP phenomenon. 

The precursors of Young-TOs thus have to be lumi- 
nous nuclei with substantial amounts of star-formation and 
possibly a bright AGN too. Another clue is that these pre- 
cursors must be found in the local-universe, since t < 10 9 
yr corresponds to z < 0.1 for any reasonable cosmology. 
Two plausible contenders for the progenitors of Young-TOs 
are starburst nuclei and starburst + Seyfert 2 composites 
like Mrk 477, Mrk 1210 and others (Heckman etal. 1997; 
Storchi-Bergmann, Cid Fernandes & Schmitt 1998; Gonzalez 
Delgado, Heckman & Leitherer 2001). Given the tendency of 
TOs to have later Hubble types than LINERs and Seyferts 
(Ho, Filippenko & Sargent 2003), it seems more attractive 
to link Young-TOs with starburst nuclei. However, the sub- 
stantial overlap in morphological properties between TOs 
and both starburst and AGN hosts, coupled to indications 
that starburst + Seyfert 2 composites have rather late type 
morphologies for AGN (Storchi-Bergmann et al. 2001) pre- 
vents us from drawing a firm conclusion at this stage. 

Similar arguments can be used to sketch the future evo- 
lution of Young-TOs. As the CYP fades, it will eventually 
cross the xy+i = 1/3 threshold below which we would not 
identify it anymore and the system would be classified as 
Old ( H4.2.H . For instance, starting from a current value of 
xy+i = 2/3, and assuming the old populations does not 
change much, the CYP would cross the xy+i ~ 1/3 line 



after it fades by a factor of 4. For an assumed age of 1 
Gyr, this would take ~ 2 Gyr to happen. In other words, 
the stellar populations of Young-TOs will become indistin- 
guishable from those of Old-LLAGN in a few Gyrs. Though 
it is tempting to link Young to Old-TOs because of their 
identical emission line properties, as pointed out in Paper II 
we cannot rule out the possibility that [OI]/Ha increases as 
the CYP fades, which would turn a Young-TO into an Old- 
LINER. Note also that for either of these two evolutionary 
connections to work Young-TOs must some how get rid of 
their excess dust ( H4.2.2I in a few Gyr, either by converting 
it to new stars or blowing it away. 

Although much work remains to be done, these gen- 
eral considerations illustrate how the careful dissection of 
stellar populations properties can provide new and impor- 
tant pieces in the quest to solve the puzzle of active galactic 
nuclei. The evolutionary scenarios sketched above will be 
examined more closely in forthcoming communications. 



5 CONCLUSIONS 

In this third paper in our series dedicated to the stellar popu- 
lations of LLAGN, we have investigated the radial variations 
of stellar populations properties in a sample of 42 LINERs 
and TOs plus 5 non-active galaxies. The analysis was based 
on high quality 3500-5500 A long-slit spectra covering an- 
gular regions of at least ~ 10" in diameter with a resolution 
of ~ 1" (corresponding to ~ 100 pc). 

The main result of Papers I and II was the identifi- 
cation of a class of objects which stand apart from other 
LLAGN in having a strong 10 8-9 yr population. In terms of 
emission lines nearly all of these nuclei have weak [OI]/Ha, 
hence their denomination as "Young-TOs". Here we have 
shown that Young-TOs are also distinct from other LLAGN 
in terms of the way stellar populations and dust are spatially 
distributed. This general conclusion was reached through 
two distinct and complementary ways. 
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First, radial profiles of absorption line equivalent 
widths, continuum colours and the total flux along the slit 
were used to trace the spatial distribution of stellar popula- 
tions. The results of this empirical analysis can be summa- 
rized as follows. 

(i) We find that the W\ profiles are of essentially two 
types: flat and diluted. Flat profiles, which indicate spatially 
uniform stellar populations, are more common, accounting 
for ~ 60% of the sample. They occur exclusively in galaxies 
dominated by an old, bulge-like stellar population, regard- 
less of the LINER/TO emission line classification. 

(ii) Diluted profiles, on the other hand, are produced by 
a central "young" population (CYP) dominated by stars of 
10 -10 yr, whose relatively blue continuum dilutes the W\s 
of metal lines with respect to their off- nuclear values. 

(iii) Although concentrated in the nucleus, these CYPs 
are spatially extended, reaching distances of up to 400 pc 
from the nucleus. 

(iv) The relation between diluted profiles and nuclear 
stellar population is clearly expressed by the ~ one-to-one 
relation between the radial dilution index 8\ and the nuclear 
W\ for the Call K line: Virtually all sources with 5k > 10% 
have W K UC < 15 A and vice-versa. This range of W K UC cor- 
responds exactly to our definition of "Young" stellar popu- 
lation, meaning populations of 1 Gyr or less. 

(v) Since these stars are found almost exclusively in ob- 
jects with [OI]/Hq < 0.25 (Papers I and II), it follows that 
stellar population gradients are typical of Young-TOs. The 
fact that these stars are located in their central regions, and 
not spread over the whole galaxy, reinforces the suggestion 
that they are somehow connected to the ionization of the 
nuclear gas. 



nucleus would also be swamped by these young starbursts, 
unless it too was brighter in the past. 

This investigation has unveiled several interesting con- 
nections between stellar population, emission line proper- 
ties, spatial distribution and extinction, paving the road to a 
better understanding of the physics of low luminosity AGN. 
Future papers in this series will explore these and other con- 
nections in further detail. 
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Second, a more detailed analysis of stellar population gra- 
dients was achieved by means of a decomposition of each 
spectra in terms of templates representative of very young 
(< 10 7 yr), intermediate age (10 8 ~ 9 yr) and old (10 10 yr) 
stellar populations. This analysis shows that: 

(vi) The CYPs in Young-TOs account for at least ~ 1/3 
of the total flux at 4020 A. We confirm the finding of Papers 
I and II that these populations are dominated by 10 8 -10 9 
yr stars. Young starbursts, even when present, make a small 
contribution to the optical light. 

(vii) Yet another property which distinguishes Young- 
TOs from other members of the LLAGN family is dust con- 
tent. Young-TOs are ~ 3 times more extincted than Old- 
LINERs and Old-TOs. This finding is confirmed using the 
HFS97 measurements of the Ha/H/3 ratio. 

(viii) Dust tends to be concentrated towards the nucleus, 
although asymmetric extinction profiles are also common. 

(ix) The radial flux distribution of CYPs have HWHM 
radii of ~ 100 pc or less. While their core is at best partly 
resolved in our data, their outer regions are clearly resolved. 

(x) The 4020 A luminosities of the CYPs are within an 
order of magnitude of 10 4 ' 3 L Q A -1 , implying B-band abso- 
lute magnitudes of ~ —15 and masses of order ~ 10 7 -10 8 
Mq. This population was 10-100 times more luminous in 
their formation epoch, at which time young massive stars 
would have completely outshone the bulge light. The active 
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